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ABSTRACT 

£N| ■ We present a closed framework of AdS/CFT with finite t/(l)s-charge chemical 

potential. We show how the gauge-invariant identification of the chemical potential 
with the bulk gauge field emerges from the standard AdS/CFT dictionary. Physical 
importance and necessity of the Minkowski embeddings within the present framework 
is also shown numerically in the D3-D7 systems. We point out that the D3-D7 model 



with only the black-hole embeddings does not have the low-temperature and low- 
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chemical-potential region in the grand-canonical ensemble, hence it is incomplete. A 
physical interpretation that explains these numerical results is also proposed. 
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I. INTRODUCTION 

AdS/CFT correspondence [1] is a useful framework to analyze strongly coupled Yang- 
Mills (YM) theories, and its application to quark-hadron physics is one of the important 
subjects from the phenomeno logical point of view. Since macroscopic properties of quark- 
hadron systems, based on thermodynamics and hydrodynamics, are as important as their 
microscopic nature, it is quite significant to establish holographic descriptions of thermo- 
dynamics and hydrodynamics of the YM theories. Finite-temperature AdS/CFT has been 
initiated in Ref. [2]. Hydrodynamic quantities of static YM-theory plasma have been com- 
puted in AdS/CFT (see, for example, reviews [3] and the references therein), holographic 
description of hydrodynamics of time-dependent YM-theory fluid has been investigated in 
Refs. [4]. 

For a complete description of the thermodynamic and hydrodynamic properties, we need 
to introduce chemical potentials of the conserved charges to the framework. Since the lattice 
gauge theory has a technical difficulty in introduction of finite baryon chemical potential, 
it is quite significant to establish a holographic description of baryon chemical potential. A 
holographic description of R-charge chemical potential has been proposed in Refs. [5, 6] and 
isospin chemical potential in AdS/CFT has also been studied in Refs. [7]. However attempts 
to introduce a baryon chemical potential to AdS/CFT have been just started recently [8-11], 
and there still exists a point which is under debate [10, 11]. The issue under the question 
is existence of so-called Minkowski-embedding phase in the D3-D7 systems [11]. For recent 
studies on finite baryon density systems in holographic frameworks, see Refs. [12, 13]. 

In this paper, we present a closed framework of AdS/CFT with finite f/(l)s-charge 
chemical potential [28]. Although our approach has wide overlap with what have been 
discussed in Refs. [8-11], the following points will be clarified in this paper: 

• A standard dictionary of AdS/CFT implies that the chemical potential is given by the 
boundary value of the zeroth component of the bulk ?7(l)-gauge field [8, 9]. However, 
this identification is not manifestly gauge invariant. A manifestly gauge-invariant 
identification of the chemical potential and the ?7(l)-gauge field has been proposed in 
Refs. [10, 11]. We clarify how the gauge- invariant formulation emerges starting from 
the standard AdS/CFT dictionary. 
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• The Minkowski embeddings at finite baryon density is claimed to be unphysical in 
Ref. [11]. However, we will show their necessity and physical significance within 
the context of our framework. It will be shown that the model with only the black- 
hole embeddings proposed in Ref. [11] has a serious problem: the model lacks the 
low-temperature and the low-chemical-potential region of the parameter space in the 
grand-canonical ensemble. We call this "incomplete-ness problem" in this paper. 

• The incomplete-ness of the model with only the black-hole embeddings can also be 
seen in terms of thermodynamic instability in the canonical ensemble. There is a 
parameter region where thermodynamically stable black-hole embeddings do not exist 
in the canonical ensemble. The Minkowski embeddings provide a stable final state in 
that case. 

• We will present a possible physical picture which clarify the difference between the 
model of Ref. [11] and that of this paper. We will also propose an idea which may 
remedy the incomplete-ness problem in the model of Ref. [11]. 

Although we will work out along the D3-D7 systems, the formalism of the C/(l) B -charge 
chemical potential given in this paper is applicable to general setups of AdS/CFT with 
flavor branes. 

The organization of the paper is as follows. In section II, we present a basic setup and 
notations in our framework. A closed formulation of AdS/CFT with finite £7(l)#-charge 
chemical potential is given in section III. The numerical results that support consistency of 
our framework is presented in section IV. The incomplete-ness of the model proposed in Ref. 
[11] is also pointed out there. In section V, we re-examine the importance of the Minkowski 
embeddings from the viewpoint of thermodynamic stability in the canonical ensemble. In 
the discussion section, we propose a possible physical interpretation of our framework. We 
also discuss how the discrepancy between Ref. [10] and Ref. [11] can be interpreted. A 
possible improvement to remedy the incomplete-ness problem of the model of Ref. [11] is 
also discussed. 
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II. BASIC SETUP 

The system we deal with is so-called D3-D7 system which corresponds to a Af = 2 SYM 
theory with flavor quarks [14]. We assume that the number of the flavors Nf is small enough 
comparing to the number of the colors N c , and we employ the probe approximation where 
the back reaction from the D7-brane to the bulk geometry is neglected. 



A. Notations 



In order to clarify the notations, we briefly review the basics of the D3-D7 system. Our 
bulk metric in the Euclidean signature is given by 

ds 2 = g (f(U)dt 2 + dx 2 ) + R 2 (j*^ + d^ , (1) 

where R 4 = 2\l A s and f(U) = 1 - (U /U) A . Here, A = g\ M N c is the 't Hooft coupling of the 
YM theory. There is a horizon at U = Uq, and the Hawking temperature is given by 

irR 2 V2\ttP s ' 1 ) 

Introducing a dimensionless coordinate £ defined by d^ 2 /^ 2 = dll 2 / (U 2 f (U)) , the bulk 
geometry becomes 

U 2 R 2 
ds 2 = — (fdt 2 + dx 2 ) + —dsl with (3) 

dsl = di 2 + i 2 dn\ = dp 2 + p 2 dttl + dy 2 + y 2 d V 2 , 

where £ 2 = y 2 + p 2 and p is the radius of the 3-sphere. U and £ are related by U 2 /Uq = 

W 2 + l/e) and/ = (l-£ 4 ) 2 /(l+e 4 ) 2 - 

Now, the induced metric on the D7-brane is given by 

ds 2 m = ^ (fdt 2 + dx 2 ) + ^ ((1 + y 2 )d P 2 + p 2 dQj) , (4) 



where y = d p y(p). 
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The DBI action of the D7-branes is 

Sdbi = N f fi 7 J dtd 3 xdpdtt 3 (g xx gnn)% \J {9ttg PP - (27va'F pt ) 2 ), (5) 

where p, 7 = [(27r) 6 Zf X/N^ 1 is the D7-brane tension and we have assumed that the non-zero 
field strength is only F pt = d p A — doA p . Notice that the Chern-Simons term vanishes within 
such a setup. Introducing u;±(£) = 1 ± £~ 4 , the DBI action is written compactly as 




Here, F pt = F pt /m T with m T = \\f\T and 

r 7 = N ff i 7 n 3 U*V 3 /4 = N f N c T 4 XV 3 /32, (7) 

where V 3 is the volume of the space where the YM theory lives. We set V 3 = 1 in this paper 
and all the extensive quantities should be understood as those per unit volume. 



B. Brief sketch of the D7-brane embeddings 

There are two big categories of the D7-brane solutions (embeddings) in the present setup: 
Minkowski embeddings and black-hole embeddings. [29] The Minkowski embeddings are the 
D7-brane configurations which do not touch the black-hole horizon, and the black-hole em- 
beddings are those connected to the black-hole horizon. 

Let us define the two foregoing embeddings by using the bulk coordinate. y(p) on the D7- 
brane worldvolume starts at a particular value at the boundary (which we define y(oo) = L) 
and one finds that it decreases monotonically as p decreases. The minimum value of y (which 
we define yo) is realized at the minimum p (which is denoted by p m in)- Since the location of 
the horizon is where y 2 + p 2 = 1, the brane touches the horizon at p = p m in > if yo < 1. 
This defines the black-hole embeddings. The Minkowski embeddings are those with y > 1 
and p min = 0. A further detailed classification of the brane embeddings is given in appendix 
A. 
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In fact, L, which is the asymptotic value of y is a parameter of the theory which describes 
the current quark mass m q . The relationship between the mass and L is given by 

m q = Uf\TL. (8) 

Since we usually fix the current quark mass and the 't Hooft coupling in the analysis, L~ L 
is equivalent to the temperature up to a proportional constant. 



III. [/(I)b-CHARGE chemical potential in ads/cft 

A. Standard AdS/CFT dictionary 

Let us consider a quark-current operator th = ip^ip in the YM theory [30]. Here, ip 
is the quark field, and we define the quark- number density Q as if;) [31]. Adding a term 
j^O^ to the YM-theory Lagrangian corresponds to switching on a one-form bulk field <3> M 
that has the following asymptotic behaviour: 

*»M = P0E + "^1 + ---. (»> 

where p — 1 for the one-form field and a is a constant which will be determined later. A is 
the conformal dimension of the operator that is 3 in this case. 

In the YM-theory side, C/(l) B -charge chemical potential /i_b in the grand-canonical en- 
semble is introduced by adding an operator joO° = [iqip^ip to the Euclidean YM Lagrangian. 
Here, jo = f^q — Hb/N c is the quark-charge chemical potential and we may use fi q rather 
than hb when it is convenient. What we need to do first in the gravity dual is to find an 
appropriate one-form field $ M . 

A proposal given in Refs. [8, 9] is that <E> M is the U(l) gauge field on the flavor D- 
branes. U(1) B symmetry is the [/(l)-diagonal part of the global flavor symmetry in the YM 
side. An important fact is that the global flavor symmetry is promoted to a gauge symmetry 
on the flavor branes in the gravity dual. For example, the global U(Nf) symmetry in the 
M = 2 SYM theory is realized as U(Nf) gauge symmetry on the flavor D7-branes. Hence 
the U(1)b symmetry corresponds to the U(l) gauge symmetry on the flavor brane, and the 
U(1) B charge is identified as the U(l) "electric" charge on the brane in terms of the gravity 
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dual. Then it is quite natural to infer that the chemical potential conjugate to the U(1) B 
charge is given by the boundary value of A that is the conjugate field to the "electric" 
charge: 



Here, we have normalized the unit electric charge so that it has the unit quark-number 
charge (that is 1/N C times the unit baryon- number charge) rather than the unit baryon- 
number charge. This normalization comes from the fact that a quark is represented by a 
single fundamental string, and the end point of the fundamental string carries unit electric 
charge with respect to A . 

Now, a problem comes up at (10); Aq is not a gauge invariant quantity in the bulk theory 
while fi q and V) in the right-hand side are physical quantities. We will provide a gauge- 
invariant definition of the quark/baryon-charge chemical potential in the next subsection. 



In this section, we construct a gauge-invariant formulation of the £/(l)£-charge chemical 
potential in AdS/CFT. A gauge-invariant definition of the chemical potential has already 
been introduced in Refs. [10, 11]. However, the aim of this section is to present a mathe- 
matically consistent framework by showing how the manifestly gauge-invariant formulation 
is obtained starting with the standard AdS/CFT dictionary (10). 



In the grand-canonical ensemble, the grand potential can be evaluated from an on-shell 
Euclidean bulk Lagrangian obtained by fixing the chemical potential at a given value. Let 
us begin with a bulk Lagrangian given by 



where we have explicitly introduced a Lagrange multiplier, denoted by Q, that defines the 
chemical potential. We omit the Lagrangian of the bulk gravity, since it is independent 



Ao(p,x) = n q + a(V>V)— H • 



(10) 



B. Gauge-invariant formulation of chemical potential 



1. Grand potential in the gravity- dual picture 




(11) 
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of the C/(l)fi-charge under the probe approximation. (11) can be re-written in a way that 
makes its AdS/CFT antecedents explicit: 



The last term is the usual bulk-boundary coupling in AdS/CFT which directly allows us 
to interpret Q as the conserved U(1)b charge of the field theory. The first two terms then 
may be viewed as the Lagrangian of the D7-branes with a charged source at p m in- We note 
that in the present background (with the only gauge field present being F pt ), the D7-brane 
Lagrangian does not have a coupling of the form QA ; thus latter term must be due to some 
other dynamical object. We will discuss the nature of the source term later. 

The third term is understood in the following way. In field theory, the grand-canonical 
partition function is defined as 



where Q is the charge and denotes the fields whose Lagrangian is L. The third term 
in the bulk action then corresponds to the field theoretic term nQ which is present in the 
grand-canonical ensemble. 

In the grand-canonical ensemble, the amount of the charge fluctuates and its expectation 
value may be determined (approximately) by the saddle point value of the partition function. 
In bulk terms, this is equivalent to determining the value of Q by minimizing the bulk action 
(11) at the given chemical potential. Thus, the saddle-point value (the on-shell value) of the 
bulk action is identified with the grand-potential of the dual field theory. 

Note that, in the usual AdS/CFT correspondence, field theory quantities are defined by 
boundary values of bulk fields. Our definition of the chemical potential in (11) is however 



To understand how this can be the case, it is worthwhile recalling how the grand potential 
in the R-charged system has been computed in terms of gravity dual. In Ref . [5] , the grand 
potential is identified with the on-shell value of the Einstein-Maxwell- Anti-deSitter effective 
Lagrangian in the bulk computed by solving the equations of motion with the boundary 



Lduai = L D7 (F pt ) + QA (p min ) +Q/i- QA (oo). 



(12) 





(13) 
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value of the time component of the vector potential kept fixed. This matches the picture 
that the control parameter of the grand-canonical ensemble, the R-charge chemical potential 
in this case, is given by the boundary value of the potential. Let us start with a slightly 
different bulk action by following this spirit: 



In this action, we have a source of strength Q at p m in, while the chemical potential has 
been defined to be the boundary value of the bulk gauge field using a Lagrange multiplier A; 
our AdS/CFT dictionary is (10) at this stage. However, demanding gauge invariance forces 
Q — A. Then, extremizing over A leads us to the definition of chemical potential as in (13). 

On the other hand, if we first extremize over A, then it is easily seen that the chemical 
potential is obtained as the asymptotic value of A , which in turn is obtained by solving the 
equations of motion of the gauge field in the presence of a charged source Q at p m i n . A = Q 
is then obtained as a consequence of flux conservation. 

The difference, of course, is that it is only in the first case that we have off-shell gauge 
invariance. In the grand-canonical ensemble, since the charge Q fluctuates, the off-shell 
gauge invariance is required. Thus we must first impose Gauss' law and then extremize 
over the Lagrange multiplier. This explains how, starting from the standard AdS/CFT 
prescription, we arrive at the definition (13). 

Indeed, the definition of the chemical potential (13) is consistent with that in Ref. [9] 
since A (p min ) is set to be zero there [32]. One interesting observation is that the right- 
hand side of (13) is nothing but a work necessary to bring a unit charge from the boundary 
to p min against the electric field F pt along the worldvolume of the D7-brane. It would be 
interesting if we can connect this picture to some physical process in the YM theory. Now 
the dictionary (10) is modified to be 



This new recipe also matches the standard computation of the grand potential in thermal 



L 'duai = L D7 {F pt ) + QA (p min ) + A(/i - A (oo)). 



(14) 



A (p,x) - A (p min ,x) = p q + a{^ijj)— H . 



(15) 
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field theories. Let us rewrite L, 



dual 



to be 



'dual 




(16) 



where 



oc 



D7 — 



dpC D7 (F pt ) + Q {A (oo) - A (p min )} 



(17) 



Pmin 



The grand potential is the on-shell value of L' D7 — p q Q obtained by solving the equations of 
motion of the fields and Q. This is nothing but the standard computational recipe of the 
grand potential if we regard L' D7 corresponds to the effective Lagrangian of the field theory. 
Indeed, L' D7 is more consistent than L D7 as a D7-brane Lagrangian since the conservation 
of the electric flux is manifest. 

In summary, the grand potential is determined by the following steps in the gravity dual: 

1. We start with L dua i. The chemical potential \i q is our input parameter. 

2. We solve the equations of motion. The boundary condition for A is given by our 
input ji q through \i q = A (oo) — A (p min ) that comes from equation of motion of Q. 

3. The on-shell value of Lduai (see appendix C for the renormalization) gives the grand 
potential (density) VL of the flavor part. The on-shell value of Q is the thermal expec- 
tation value of the quark-charge density. 

We employ this gauge-invariant formulation of the chemical potential in the present paper, 
and we choose A p = gauge in the remaining part. 
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2. Connection to thermodynamic relation 



A connection between the equations of motion and a thermodynamic relation can also 
be seen in the following way. The equations of motion of A from (11) are 



d ( dCo7 



dp V dA 
dC D7 



A 
dC D7 



dA 



o 



p=oo 



P — Pmin 



o, 

-Q, 
-Q. 



(18) 
(19) 
(20) 



Solving the differential equation (18) with boundary conditions (19) and (20), we obtain 



dC 



D7 



dA 



-Q, 



(21) 



which suggests 5 Cm = — QSA . Since Q does not depend on p, we have 

/OO PCX) 
dpC D7 = -Q 5 / dpA , 
, min J Pmin 

which is equivalent to the thermodynamic relation 



(22) 



on 



dp q 



-Q. 



(23) 



Here we have used the definition of the chemical potential (13) and the fact that the on-shell 
value of Ljuai is given by that of f°° dpCm- It is interesting that the Gauss-law constraint 
(21) in the "electro-magnetism" on the flavor brane corresponds to the thermodynamic 
relation (23) in the YM theory. 

The coefficient a in (9) can also be determined by the boundary condition (19). Substi- 
tuting the asymptotic expansions (15) and y = L + 0(p~ 2 ), 



-Q 



dC 



D7 



A 



p=oo 



--T 2 N c N f a(^). 



(24) 



To make the foregoing expression to be consistent with our identification Q = (ip^ip), a is 
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determined to be 

a= wh;- (25) 

C. Legendre transformation to the canonical ensemble 

The Helmholtz free energy F per unit volume is obtained by the Legendre transformation 
F = Q + p q Q. This is nothing but a procedure to obtain an on-shell "Hamiltonian" (if we 
regard p as a "time" formally) in terms of the D7-brane dynamics: 



/oo 
dp 
, min 



oA 



(26) 

on— shell 



Here in the right-hand side, — Q is interpreted as the "conjugate momentum" of A by (21), 
and all the quantities in (26) are the on-shell values. 

Then, we can make the following recipe in the canonical ensemble. 

1. We define the Hamiltonian 

H(Q) = C D7 - -f^A , (27) 

where Ti and Cr,7 are not necessarily at on-shell, but the Gauss-law constraint (21) is 
imposed to eliminate the electric field in favor of Q. Now our input parameter is Q 
but not p q . 

2. We solve the equations of motion to find the on-shell value of H which gives F. 

3. The thermal expectation value of p„ is given by f °° dpF pt , where F pt is obtained from 
the Gauss-law constraint (21). 

The Hamiltonian for our system is explicitly given by 

H = V(y,p)^/TT¥, withV(y,p) = rJ^(^Q^ + colp^. (28) 
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The equation of motion after eliminating the gauge field is written explicitly by 



y 



+ 



d log V . d log V 



= 0, 



(29) 



1 + y 2 



dp dy 



and the Gauss-law constraint in the explicit form is 



F p t = Q 




(30) 



Notice that V(y, p) can be regarded as an effective tension of the D7-branes which is affected 
by the presence of the gauge field. In particular, the effective tension becomes larger if the 
charge increase, and the brane is attracted to the black-hole more strongly. Therefore if we 
keep the current quark mass fixed, the minimum "height" yo of the D7-brane decreases as 



both V(y, p) and (30) are written in terms of Q. We may use Q as a parameter of the system 
in the canonical ensemble in the later analysis if it is convenient. 

Before closing this section, we would like to make a few comments. The Gauss-law 
constraint in the form of (30) gives a relationship between p q and Q: an equation of state. 
Equation (30) itself can be used both in the grand-canonical ensemble and the canonical 
ensemble. We should be careful only in the interpretation of the variables: in the (grand-) 
canonical ensemble, p q (Q) obtained from (30) is the thermal expectation value while Q (p, q ) 
is the control parameter of the theory. 

Another comment is about the framework presented in Ref. [10]. There is a "quick" 
computational method based on our formalism. Since the on-shell value of L^uai is the same 
as that of L_D7, the correct grand potential is also obtained by starting with Ld7 alone, and 
solving the equations of motion by imposing the Gauss-law constraint ^j 1 - = —Q by hand. 
This is nothing but the computational method employed in Ref. [10]. This means that Ref. 
[10] produces the results of the present framework. 

The last comment is about the consistency of the present formalism. We would like to 
stress that the charged source whose necessity is pointed out in Ref. [11] has already been 
introduced at the stage of (12). Thus, there is no need for introducing any further sources, 



we increase the charge density. It may be useful to define Q = 



16 



Q since 



^/XN c NfT 3 
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and all equations are correctly satisfied. 

IV. CONSISTENCY IN NUMERICAL RESULTS 

Now, we would like to bring our attention to the claim presented in Ref. [11]. The authors 
of Ref. [11] pointed out the necessity of the charged source at p = p m in, an d they provided 
it by adding U(1)b- charge carrying objects to the system. Their idea is to put fundamental 
strings (Fl's) between the D7-branes and the black-hole horizon, that are interpreted as 
quarks. Then, they found that the Minkowski embeddings in Ref. [10] are unstable due to 
the tension of the Fl's and there is no way to keep the D7-branes off the horizon. This is why 
the Minkowski embeddings at finite baryon-number density are concluded to be unphysical 
in Ref. [11]. 

However, one should notice that this is not what we are doing in the present paper. The 
necessary charged source term has been introduced as an external source to the D7-brane 
DBI theory at (12) and we have not added any corresponding Nambu-Goto action of the 
Fl there. This means the system we are dealing with is something different from that in 
Ref. [11]. Since we have not introduced the additional Fl's, we expect that in the present 
framework, the Minkowski embeddings are physical as well. 

In this section, we will show that it is indeed the case. We will present numerical results 
[33] which indicate consistency of the present formalism and the necessity of the Minkowski 
embeddings. A physical interpretation of our setup will be proposed in section VI. 

A. p-Q diagram and Maxwell construction 

In order to show a consistency of our framework, we will employ p-Q diagrams where 
the relationship between p q and Q obtained from (30) is drawn. (We may use p as the 
meaning of p q .) Let us present basic explanations on the p-Q diagram in this subsection as 
preparation. 

An example of the p-Q diagram is given at Fig. 1(a) (or at Fig. 1(c)) [10]. The D7-brane 
solutions that belong to the Minkowski embeddings start from A (which is the origin of 
the plane) and go through B and C until D where the line meets the vertical axes. The 
black-hole embeddings exactly start from D where the Minkowski embeddings terminate, 
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and go through E, F, G, H, I and extend to the large-Q and large-// region. 

It is worthwhile checking a consistency of the fx-Q diagram. Fig. 1(b) and Fig. 1(d) 
show the relationship between the Helmholtz free energy density and Q, and that between 
the grand potential density and //, respectively. Although the thermodynamic potentials 
obtained from all the possible solutions are indicated, what we should take is the line which 
has the minimum value. Then, we find first-order phase transitions. Let us see Fig. 1(d), 
for example. The phase transition is a jump between a Minkowski embedding and a black- 
hole embedding at the critical chemical potential fj,i. Since the grand potentials in both 
embeddings have the same values at /ii, the integral f Qd/i below /ii and that above /xi 
have to be same because the integrals compute — VI by virtue of ^ = —Q. This is the 
Maxwell construction. One can see that the Maxwell construction works in a non-trivial 
way in good accuracy in Fig. 1(c). We have checked numerically that the area of the shaded 
regions marked "+" agrees with the area of the regions with "— ", due to a non-trivial 
collaboration of the Minkowski and the black-hole embeddings. We can also see that the 
Maxwell construction works in the canonical ensemble in Fig. 1(a), again by virtue of the 
collaboration of the two types of the embeddings. 



B. Parameter space in the grand-canonical ensemble 

Now, we are ready to make some important comments based on the numerical results. 
Fig. 2 shows that how the \x-Q diagram is deformed if we make the temperature lower. The 
temperature goes down as we move from Fig. 2(a) to Fig. 2(d). An important feature is 
that the cusp H located at the origin on Fig. 2(a) "goes up" (Fig. 2(b)) and disappears 
(Fig. 2(c) and Fig. 2(d)) along the cooling process. The cusp on Fig. 2(b) gives the 
minimum value of the chemical potential within the black-hole embeddings. This means 
that the low-chemical-potential region disappears from the parameter space of the theory at 
the sufficiently low temperature if we abandon the Minkowski embeddings. 

We can also interpret the behaviour of the diagram in a different way. Suppose that 
we are in the grand-canonical ensemble and we examine a process in which we vary the 
temperature with maintaining the chemical potential. Then, we encounter a problem that 
there is no low-temperature region in our parameter space at sufficiently small chemical 
potential. 
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1=0.764117 




HO. 00250. 0050. 0075 0.01 0.012 50. 015 



(c) (d) 

FIG. 1: (a) An example of [i-Q diagram. The shaded regions indicate the Maxwell construction 
in the canonical ensemble, (b) Determination of the critical densities from F. We observe two 
phase transitions in this case, (c) The same [i-Q diagram can be used in the grand-canonical 
ensemble. The shaded regions indicate the Maxwell construction in the grand-canonical ensemble, 
(d) Determination of the critical value of [i from Q. 

The lack of the low-temperature region can be seen more clearly in Fig. 3, where the 
D7-brane solutions on the L-y plane are given. The Minkowski embeddings (j/o > 1) 
and the black-hole embeddings (yo < 1) are connected at y = 1 there. Notice that all 
the temperature region is covered by using the two types of embeddings. However, if we 
abandon the Minkowski embeddings, the theory cannot cover the low-temperature region. 
(Recall that L ~ l/T.) 

We can conclude from the above that if we abandon the Minwkowski embeddings, the 
low-temperature and the low-chemical-potential region in the grand-canonical ensemble does 
not exist in the formalism; we have no way to introduce the flavor degree of freedom in that 
region. This is an incomplete-ness of the formalism. On the other hand, if we include the 
Minkowski embeddings together with the black-hole embeddings, all the parameter region 
of the theory is covered in harmony of the two types of embeddings. We can also see that 
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FIG. 2: [i-Q diagrams at various temperature. 
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FIG. 3: The structure of the D7-brane solutions at /x = 0.01 fixed. Solid line: black-hole embed- 
dings, dashed line: Minkowski embeddings. 



the two embeddings are connected into a single family of the solutions on Fig. 1, Fig. 2 and 
Fig. 3. We present some numerical results in the grand-canonical ensemble in appendix B 
to make the above points vivid. 
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Incidentally, we can understand why we did not encounter the above "incomplete-ness" 
problem manifestly in Ref. [11]. Let us look at Fig. 2 again. We notice that the entire 
region of Q can always be covered even if we use only the black-hole branch. This means 
that if we examine the process in which we vary the temperature with fixing Q, we have 
always at least one solution inside the black-hole branch at any temperature; all the region 
on the Q-T plane can be covered even if we use only the black-hole embeddings. Namely, we 
do not have a missing parameter region in the canonical ensemble. The reason why we do 
not see the incomplete-ness in Ref. [11] in a manifest way is that the analysis is given in the 
canonical ensemble there. However, our claim is that the canonical ensemble with only the 
black-hole embeddings is transformed to the incomplete formalism of the grand-canonical 
ensemble through the Legendre transformation: the formalism lacks something necessary 
even in the canonical ensemble. In the next section, we will see on this point and we will 
find the Minkowski embeddings play an important role in the canonical ensemble as well. 



V. THERMODYNAMIC STABILITY IN THE CANONICAL ENSEMBLE 

We will further look into the importance of the Minkowski embeddings by investigating 
the thermodynamic properties of the system in this section. A key issue here is thermody- 
namic stability of the system. We will examine the thermodynamic stability based on the 
following two conditions in the canonical ensemble: 



1 

L - dQ 



> 0: 

T 



If this condition does not hold, the system can lower its free energy by separating into 
two phases with densities Ql < Q < Qh [H]; the baryon-charge density tends to be 
inhomogeneous. We call this instability "number-density instability" in this paper. 
Existence of this instability in the D3-D7 systems has been originally discovered in 
Ref. [10]. 

Z - 8T \Q - U - 

S denotes the entropy density of the flavor part. This condition is equivalent to the 
positivity of the specific heat (in the canonical ensemble), hence it has to hold. 

The foregoing conditions have also been examined with the black-hole embeddings in Ref. 
[11] based on the canonical ensemble. The number-density instability has been found, while 
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the instability based on the condition 2 has not been found there. A detailed study on 
thermodynamic instability at zero charge density is found in Ref. [16]. 

In the present setup, we will also find the number-density instability. However, we will 
see that the Minkowski embeddings provide a stable final state. We will also comment on 
the instability based on the condition 2. 



A. Number-density instability and phase structure 

Let us go back to the \i-Q diagram, Fig. 1(a), and examine the number- density stability 
in the canonical ensemble. There is a region where < between points D and G [10]. 
(We call this region "black- hole B branch." See for the details, appendix A.) The region 
between E and F has the lowest free energy if we assume the charge density is uniformly 
distributed. However, the free energy may be lowered by allowing the system to be a mixture 
of two phases of different densities, as it is discussed in Ref. [11]. 

Let us remind us of the necessary conditions to form a stable final state as a mixture of 
two types of domains of different phases. Since the charge can move from one domain to 
the other, the condition for the stability of each domain is what we may usually use in the 
context of the grand-canonical ensemble: the chemical potentials of the two domains have to 
be equal to each other in order to achieve thermal equilibrium, and the condition > 

has to hold in each domain to make the thermal equilibrium to be stable. 

In the present setup, we have two first-order phase transitions at Q = Qi (from E to B) 
and at Q = Q 2 (from F to / ) where the system jumps into the stable phases. Since the 
chemical potential at B is larger than that of at /, the final value of the chemical potential 
should be somewhere between the two, namely the points indicated by B' and /' in Fig. 
1(c). Therefore, the final state can be realized as a stable mixed phase of the Minkowski 
phase at B' and the black-hole C phase (see for the details, appendix A) at 

However, if we do not have the Minkowski embeddings as the model in Ref. [11], the 
low-density region is not bounded by any stable phase, and one of the domains which may 
form the mixed phase is still in the unstable black-hole B branch. This is a problem: there 
is no way to achieve the thermodynamic stability in this case. 

Indeed, the problem is seen as absence of phase transition in the grand-canonical ensem- 
ble. Points B' and I' in Fig. 1(c) which forms the stable mixed phase in the canonical 
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ensemble are nothing but the points we have a phase transition from the Minkowski to 
the black-hole embeddings in the grand-canonical ensemble. If we remove the Minkowski 
embeddings, we have no phase transition in the grand-canonical ensemble. 

It would be useful to present the phase structure of the system both in the canonical 
and the grand-canonical ensemble to overview the above mentioned properties. The phase 
diagrams in the canonical ensemble are shown in Fig. 4. (See also Ref. [10].) As it has been 




FIG. 4: (a) The phase diagram in the canonical ensemble. The shaded region is where the number- 
density instability presents, (b) The region of two phase transitions exist is zoomed. 



reported in Ref. [10], there are two phase transitions in the small region of 

Q[ = 4.9 x 10~ 4 < Q < = 3.9 x 10~ 3 . (31) 

The phase transitions are of the first order except for the endpoint of the upper line at Q = Q* 
where a second-order phase transition takes place. A part of the phase diagram is zoomed 
in Fig. 4(b). The detailed analysis shows that the three different phase boundaries meet at 
Q = Q* with different slopes. Here, the low-temperature phase is the Minkowski phase, the 
phase between the two transitions is the black-hole B phase and the high-temperature region 
is the black-hole C phase. (See appendix A, for the detailed classification of the phase.) The 
region where the number- density instability exists is indicated as a shaded region in Fig. 4. 
We can see that the unstable region is surrounded by the stable regions. This means that 
the system can reach a final stable mixed phase in any case. If we remove the Minkowski 
embeddings, the Minkowski phase disappears. The lower phase boundary disappears and 
we have only one phase transition (this is essentially what is reported in Ref. [11]). In this 
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case the unstable region reaches the T = axes and it is not bounded by the stable region 
[11]. 

The phase structure in the grand-canonical ensemble is presented in Fig. 5. A significant 



. 55 




FIG. 5: Phase diagram in the grand-canonical ensemble 

difference from the canonical-ensemble case is that we have only one phase transition, as 
one may expect from the discussion based on Fig. 1(c) and Fig. 1(d). The low-temperature 
phase is the Minkowski phase while the high-temperature phase is the black-hole phase. 
The phase transition is of the first order and the value of Q jumps at the transition point; 
the unstable region we found in the canonical ensemble is "skipped" by the jump at the 
transition (see Fig. 1(c) again). We did not find any number-density instability within 
our analysis in the grand-canonical ensemble. If we remove the Minkowski embeddings, the 
Minkowski phase disappears and the system has no phase transition. Furthermore, the low- 
temperature low-chemical-potential region will not be covered by any embeddings. We can 
see the cooperation of the Minkowski and the black-hole embeddings in the grand-canonical 
ensemble vividly in appendix B. 



B. The specific heat and the third law of thermodynamics 

Here, we would like to comment on the temperature dependence of the entropy density 
S. We may use 

ST = -AF + cm q + 3Qfi, (32) 
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rather than S 



dF{T,Q) 
8T 



where F is the Helmholtz free energy density computed only 



from the DBI theory of the flavor D7-branes and the quark condensate. (32) is 

derived by using the fact that the free energy is a function of m q /T and Q/T 3 in the form 



of F = T 4 f(m q /T, Q/T 3 ). (32) may be technically more useful than S 



8F(T,Q) 



8T 



in the 



numerical analysis since we need not to differentiate the numerical data. Notice that the 
above entropy density is not that of the total system. Since we are not taking account of 
the bulk (adjoint fields) free energy, the entropy we discuss here is that of the flavor part. 

A typical, but preliminary result on the temperature dependence of the entropy density 
in the large-Q region is plotted in Fig. 6. 
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FIG. 6: (a) S as a function of T and (b) S/T 3 as a function of T, at Q = 1.0 fixed. Diamond-shaped 
plots with solid blue line: black-hole embeddings, box and star-shaped plots with dashed red line: 
Minkowski embeddings. The vertical lines indicate the position of T c . 



We find that the entropy still remains finite at T = in the black-hole phase at this 
parameter region. This means that the third law of thermodynamics will be broken if we 
abandon the Minkowski embeddings. Indeed, S/T 3 that counts the degree of freedom per 
unit volume of the system diverges at T = in the black-hole phase. However, this is not 
physically realized in our setup since the system at T = is in the Minkowski phase where 
the third law of thermodynamics seems to hold. This is another example of the role of the 
Minkowski embeddings. 

One may notice that we still have a region of J^L < in the black-hole phase even if we 
take the Minkowski embeddings into account. The specific heat of the flavor brane in this 
region is negative and the flavor part is not at a stable thermal equilibrium. However, the 
entropy density we analyzed is only that in the flavor part; a stable phase may be realized 



23 



after interaction of the flavor part and the bulk part. This implies that we may need to 
go beyond the probe approximation. (This also matches the fact that this instability is 
observed only in the large-Q region where the D7-brane's effective tension V(y,p) becomes 
large.) Because of this reason, the statement in this subsection should be understood as a 
preliminary comment. We leave detailed study on this negative-slope behaviour for future 
work. 

VI. DISCUSSION: AN INTERPRETATION OF THE PRESENT SETUP 

We have presented a gauge-invariant formulation of U (1) ^-charge chemical potential in 
AdS/CFT. The framework was constructed in a mathematically closed way. We have also 
numerically examined the properties of the system. The results indicate that it is consistent 
to employ all the possible types of the D7-brane embeddings in our framework, and the 
necessity of the Minkowski embeddings has been shown. Then, how can we reconcile the 
claim in Ref. [11] with our formalism? The purpose of this section is to present one possible 
physical interpretation of our framework, and to clarify the difference between our system 
and that of Ref. [11]. 

To this end, let us consider the Lagrangian 



The second term which is the source term does not arise from any known coupling of the 
D7-brane (in this background). The remaining terms however, can be argued for, from the 
standard AdS/CFT correspondence (as was done earlier). Thus, one is forced to regard the 
second term as coming from a charged object which intersects the D7-brane at p m in- 

The authors of Ref. [11] have asked what is the physical object carrying the U (1) ^-charge 
in the gravity-dual side. Their proposal is that this object is a fundamental string connecting 
the D7-branes and the black-hole horizon which corresponds to a single quark in the YM 
theory. Hence they added the fundamental strings (Fl's) further to the system. 

The bulk action for the system in Ref. [11] in our language is then, 



Lduai = L D7 + QA (p min ) +Q(fi- A)(oo)). 



(33) 



LrMMMT — Ldual + QV^Lng, 



(34) 
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where L NG is the Nambu-Goto Lagrangian of the Fl, and QV 3 L NG represents the contribu- 
tion of the bundle of the Fl's that are homogeneously distributed in the three-dimensional 
space with number density Q [11]. In this picture, the second term in (33) is induced by the 
end points of the Fl's on the D7-branes. Now the equations of motion of the system will 
be modified because of QV^L^g- in particular, the effect of the tension of the bundle of the 
Fl's has to be taken into account (this force balance condition may be regarded as arising 
from the surface terms in the equations of motion). 

The conclusion of Ref. [11] is that the bundle's tension is always strong enough to pull 
down the D7-branes to the horizon, and we have only black-hole embeddings if we start 
from (34). However, this is clearly different from what we are doing in this paper. Since the 
Nambu-Goto Lagrangian is absent from our action (33), there is no force-balance problem 
between the D7-branes and the Fl's, and there is no reason to abandon the Minkowski 
embeddings a priori in our framework. Then what is the physical interpretation of our 
setup? Let us look into (34) in more details to clarify the physical interpretation. 

If we start from (34) and carry out constructing the formalism (by regarding Q as a 
Legendre multiplier), the definition of the chemical potential will be modified to 



The mass of the Fl (which is defined to be L NG here) is now included in the definition of 
the chemical potential. 

Here we should notice that we could have added baryons instead of quarks to the system a 
la Ref. [11]. A baryon is represented by a D5-brane wrapped on the S 5 with N c fundamental 
strings attached [19], which is often called a baryon vertex. (See also Refs. [20].) If we add 
(a bundle of) the baryon vertices instead of the Fl's, the second term in the right-hand 
side of (33) is understood as a term induced by the Chern-Simons term of the D5-brane 
action. We should also add the "baryon- vertex Lagrangian" (D5-brane DBI Lagrangian and 
Nambu-Goto Lagrangians of the fundamental strings that constitute the baryon vertex) to 
(33) instead of the Fl's Nambu-Goto Lagrangian: 



jj q + L NG + Q 




(35) 



L 



'new 




(36) 



c 
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where L baryon is the Lagrangian of the baryon vertex. Now the definition of the chemical 
potential will be modified to be pi new = fi q + + However, this is not what 

we are doing, either. 

The point is that the definition of our chemical potential [x q is blind to what the baryon- 
charge carrying object is. Therefore, one natural interpretation of our formalism is that we 
are dealing with a system of mesons with an external {7(l)£-charged source in the canonical 
ensemble (and its conjugate system in the grand-canonical ensemble). This is why we 
have used an expression "C/(l) B -charge chemical potential" rather than "baryon chemical 
potential" in this paper. This interpretation comes form the fact that the D7-brane's DBI 
action provides an effective action of mesons, and we have inserted only the charged source 
term to it in the gravity-dual side. 

One may wonder how to imagine the system of mesons with an external ?7(1 ^-charged 
source. An intuition may be obtained by considering so-called Walecka model (or the <j-oj 
model) [21]. The Walecka model is a model of nucleon-meson systems. The baryon-number 
current couples to the cu-meson through the Yukawa coupling in that model. Therefore, the 
baryon-number density (which is charged under the U(1)b symmetry) can act as a source 
for the u° mesons there. We will not try to connect the present setup with any particular 
phenomenological model of the SYM theory in this paper. However, our baryon-number 
current can couple to our mesons in the same way as it does in the Walecka model, in 
principle. [34] 

Let us re-interpret what we have observed based on the above interpretation. First of 
all, the existence of the physical Minkowski embeddings is explained. Since we are dealing 
with a system of mesons with an external £7(l)s-charged source inserted, we need not to 
add any [/(l)s-charge carrying dynamical object to the system by hand. Therefore, the 
fundamental strings that make the Minkowski embeddings to be unstable can never come 
into the present setup by definition. 

The instability of the Minkowski embeddings discussed in Ref. [11] may also be inter- 
preted in the following way. The color non-singlet quarks have been added to the phase 
(Minkowski branch) where the physical mesons exists. The instability of the Minkowski 
embeddings with Fl's may correspond to the instability of a confinement phase with color 
non-singlet quarks. 
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The reason why the framework in Ref. [11] is incomplete may also be understood. If 
we wish to analyze the system which contains physical baryon-number charged particles, we 
also have to examine what happens if we add the baryon vertices instead of the Fl's to the 
Minkowski embeddings. This analysis has not yet been done, and what we should do at this 
stage is to postpone the conclusion about the fate of the Minkowski embeddings rather than 
to abandon them, until the analysis with baryons is accomplished. 

An argument that the charged source can be D5-branes is the following. Let us imagine 
that we introduce charged sources keeping the D7-branes fixed. We can then compare the 
energies of various candidate sources. That is to say, we compare the energy of N c Fl's 
stretching between (p, y) = (0, yo) and the horizon, and the energy of a single D5-brane 
(wrapping the S 5 ) placed at (0,yo). We can immediately see, that if yo is sufficiently large, 
D5-branes are lower in energy[18, 20]. Thus, baryons are more appropriate as charged sources 
for low enough temperature (i.e., large enough y ). (A related argument that D5-branes with 
attached fundamental strings can sit outside the horizon in the black-hole geometry for low 
enough temperatures has been demonstrated in [20]). 

Another candidate source is Dl-branes that are transverse to the D7-branes (i.e., wrap 
the <p direction of the metric (4)). It can be argued that the open strings from the Dl-branes 
to the D7-branes are fermionic [19]. Further, the Dl-branes can form a fuzzy S* 5 by a version 
of the Myers' effect [23]. 

Indeed we have a preliminary observation of what happens if we add baryon vertices to 
the Minkowski embeddings. In Fig. 7(a), the free energy densities with/without adding 
the D5-brane mass are shown there. Fig. 7(b) shows the phase diagram in the canonical 
ensemble with/without adding the D5-brane mass. The D5-brane's mass is computed based 
on the assumption that the configuration of the D5 is spherical and it is attached to the 
Minkowski embedding at (p, y) = (0, yo) without taking account of interactions between the 
D7-brane (see, for the details, appendix D). The correction to the black- hole embeddings is 
zero under this assumption since the D5-brane is on top of the horizon. This is why the line 
of the phase transition between the black-hole embeddings (between the B-branch and the C- 
branch) is not modified in Fig. 7(b). We understand that the assumption employed here may 
be invalidated if we analyze the dynamics of the full system of the D5's and the D7-branes. 
However, what we would like to comment from Fig. 7 is that the naive estimation of the 
D5-brane's energy is not too large to make the phase transition impossible; the modified free 
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energy of the Minkowski embeddings still intersects that of the black-hole embeddings. This 
suggests that we have a chance to obtain a complete setup for the finite density systems 
of dynamical baryons, by including the baryon-vertex D5-branes in an appropriate way. 
Notice that our estimation of the D5-brane mass qualitatively matches what is discussed in 
Ref . [13] ; the D5-brane mass is a decreasing function of Q since yo is. Improvement of the 
framework in this direction is certainly an important subject. 




(a) (b) 

FIG. 7: (a) F-T diagrams with/without the D5-brane correction at Q/T = 1.0. The blue solid line 
represents the black- hole embeddings, the red dashed line represents the Minkowski embeddings 
without the D5-brane mass and the purple dotted line is the Minkowski embeddings after adding 
the D5-brane mass, (b) Phase diagram in the canonical ensemble with/without the D5-brane 
correction. The blue line (the upper one) indicates the phase transition within the black-hole 
embeddings. The red line in the middle represents the phase transition between the Minkowski 
and the black-hole embeddings without taking account of the D5-brane mass. This phase boundary 
is modified to be the black line (the lowest one) after adding the D5-brane mass. 

Finally, we would like to comment on the location of the charge on the D7-branes. We 
have assumed that the location of the charge on the D7-branes is p = p m in at (12). However, 
there is a freedom in the choice of the location in principle. For example, the charge 
can be delocalized in the p-direction, or the position of the charge can be different from 
P — Pmin- If we take the picture where we have dynamical charged objects, the distribution 
of the charge in the p-direction may be determined by the dynamics of the D7-branes and 
the charge-carrying objects. However, since the charge in the this paper with the above 
interpretation is non-dynamical, the distribution should be given by hand in the present 
setup. One may wonder how to determine the location of the charge in the p-direction and 
what is its physical interpretation. In this paper, we have taken the most simplest choice 
for the charge distribution, and we leave these problem [35] open. All the points discussed 



28 



in this section have to be clarified more in further studies. 



Note added: While the manuscript of the present paper is prepared, we received related 
papers [24], [25], [26] and [27]. 
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APPENDIX A: DETAILED CLASSIFICATION OF BRANE CONFIGURATIONS 

The brane embeddings in our system have rich variety. We will present the detailed 
classification of the brane embeddings and the phases based on the canonical ensemble [36] 
[37] in this appendix. 

Minkowski embeddings 

The key feature of the Minkowski embeddings is captured by Fig. 8, where we plot yo, 
the value of y at p = 0, against the asymptotic height L = y\ p=00 - 
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FIG. 8: L vs. yo for the Minkowski embeddings for various Q. 
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In the case of zero temperature where we do not have black hole in the bulk, the equations 
of motion become very simple. One can easily show that the obtained brane configuration 
is just horizontal, namely y = yo = L, regardless of the value of Q. Let us consider more 
general setup for finite temperature cases. For large values of yo, the branes are hardly 
affected by the black hole and the brane profile is pretty much similar to the horizontal 
configuration of the T = case. This means that L decreases linearly (with the slope ~ 1) 
regardless of the value of Q as we decrease y at large y (see the large-y region in Fig. 8). 
However, there is a critical value at which L reaches a minimum. For smaller values of 
yo, the branes seem strongly "repelled" by the black hole and L increases again steeply. 

We will refer to those Minkowski embeddings in the region of yo > y$ as Minkowski A 
branch and those in the region of yo < y^ as Minkowski B branch. These two branches exist 
even when Q = 0. The A-branch always has lower free energy than the B-branch, hence the 
B-branch will not be physically realized. 

As we increase Q, y^ moves to a larger value and the corresponding minimum value of 
L also increases. This means we do not have any Minkowski embeddings in the region of 
sufficiently large Q and sufficiently high temperature, in the canonical ensemble. 

Black-hole embeddings 



The relationship between y and L for the black- hole embeddings is shown in Fig. 9. 
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FIG. 9: L vs. yo f° r the black-hole embeddings for various Q. 

As we increase yo from zero, L initially increases linearly with the common slope regardless 
of the value of Q. For very small densities Q < Q* = 8.9 x 1CT 3 (that corresponds to 
Q < Q*), L attains a local maximum (at yi ma x), and then decreases to a local minimum (at 



30 



Dimin)- This cubic behaviour is reminiscent of the Van der Waals P — V diagram. The local 
maximum of L is present even for Q = 0. As we increase Q, yi max becomes larger while yi m in 
goes smaller, and they merge at the critical value of Q = (that corresponds to Q = Q*). 
This is again analogous to the behaviour of the Van der Waals P — V curves with Q playing 
the role of temperature. 

When y becomes nearly 1, L increases sharply. In the language of configuration, the 
worldvolume starts from the horizon and grows rapidly in the y direction within a short 
range in p like the Minkowski B-type branes. 

We will classify the black-hole embeddings into three branches when we have the local 
maximum/minimum. Those in the region of y < yi max are called as black-hole C branch. 
Those in the region yi ma x < y < yimin will be referred to black-hole D branch, and those 
with yimin < y < 1 is defined as black-hole B branch. If we do not have the local maxi- 
mum/minimum, there is no further classification of the black hole embeddings. (We may 
call them as black- hole A branch.) The D-branch always has higher free energy than that 
of the black-hole B or C branch, hence it will not be realized as a physical phase. The 
black- hole B branch can have the lowest free energy among them in some case, however it 
has the number- density instability discussed in section V A. 



APPENDIX B: MORE ABOUT THE GRAND-CANONICAL ENSEMBLE 

We present some more numerical results in order to make the collaboration of the 
Minkowski and the black-hole embeddings vivid. Temperature dependences of various quan- 
tities are shown in Fig. 10. Here, the entropy density (of the flavor part) is computed by 
using 



b = ~df 



-4Q-Qfi + mqC ^ B1 ^ 



T 



where c is the quark condensate. Each diagram is represented as a connected line within the 
resolution of the numerical data that consists of the Minkowski branch (low-temperature 
region) and the black- hole branch (high-temperature region). Because of the cusp-shaped 
structure (see Fig. 10(b), for example), the physical value of the entropy density, the number 
density and the quark condensate jump at the critical temperature of the (first-order) phase 
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transition. We will not get into the details of the physical interpretation of the results in 
this paper, but one can see how the Minkowski and the black-hole embeddings cover the full 
temperature region from these examples. 



APPENDIX C: REGULARIZATION 

The on-shell value of Ld 7 (Q) is in fact divergent due to the integral over p from p min to 
infinity. This IR divergence in the bulk is understood as a UV divergence of the boundary 
theory and therefore we need to renormalize it. 

We introduce a cut-off p max . One can check that the integrand CmiQ) behaves 

^ = pL, + o( p -L). (ci) 
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The renormalization scheme we adopt is the minimal subtraction, namely we add the counter 
term 



^counter 



-TlPtnaJ^ (C2) 



to L B7 (Q). Indeed, our renormalization method is the same as that of the holographic 
renormalization given in Ref. [17]. 

APPENDIX D: BARYON MASS 

The mass of single D5-brane spherically wrapped on the S 5 is given by 

M B = ^gTt^VDs, (Dl) 

where // 5 = [(27r) 4 /gA/A^ c ] _1 is the tension of the D5-brane and V D5 = ir 3 R 5 is the volume of 
the S* 5 . Substituting the metric (1) and the constants into (Dl), we obtain [22] 

1 . at 

M b = - a U b ^J(Ub)^ (D2) 
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where U B is the location of the baryon vertex in the U direction. It is convenient to express 
M B by using the (p, y)-coordinate. Notice that we have the following relations: 

vWb) = 4^, (D3) 

Un 



U B = ^yl + y- \ (D4) 

where we have assumed that the location of the D5-brane on the p = axes is yo- Plugging 
Uq = \/2\k1 2 s T into the formulae, we obtain 

M B = ^V\TN c I l7 (D5) 

/^y^(f^). (D6) 
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The correction we added to the Helmholtz free energy of the Minkowski embeddings in 
section VI is (Q/N C )M B . 
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